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[~{ ' Summary. Except in a few cases cosmic dust can be studied in situ or in terrestrial 

, laboratories, essentially all of our information concerning the nature of cosmic dust 

depends upon its interaction with electromagnetic radiation. This chapter presents 
ON ■ the theoretical basis for describing the optical properties of dust - how it absorbs 

and scatters starlight and reradiates the absorbed energy at longer wavelengths. 
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Ql|. 1 Introduction 

O ' Dust is everywhere in the Universe: it is a ubiquitous feature of the cosmos, 

-4— > . impinging directly or indirectly on most fields of modern astronomy. It occurs 

^ ' in a wide variety of astrophysical regions, ranging from the local environment 

of the Earth to distant galaxies and quasars: from meteorites originated in 
the asteroid belt, the most pristine solar system objects - comets, and strato- 
spherically collected interplanetary dust particles (IDPs) of either cometary or 
' asteroidal origin to external galaxies (both normal and active, nearby and dis- 

£Nl , tant) and circumnuclear tori around active galactic nuclei; from circumstellar 

envelopes around evolved stars (cool red giants, AGB stars) and Wolf-Rayet 
stars, planetary nebulae, nova and supernova ejecta, and supernova remnants 
OO ' to interstellar clouds and star-forming regions; from the terrestrial zodiacal 

cloud to protoplanetary disks around young stellar objects and debris disks 
around main-sequence stars ... 

Dust plays an increasingly important role in astrophysics. It has a dramatic 
effect on the Universe by affecting the physical conditions and processes tak- 
ing place within the Universe and shaping the appearance of dusty objects 
' (e.g. cometary comae, reflection nebulae, dust disks, and galaxies) (i) as an 

absorber, scatterer, polarizer, and emitter of electromagnetic radiation; (ii) 
as a revealer of heavily obscured objects (e.g. IR sources) of which we might 
otherwise be unaware; (iii) as a driver for the mass loss of evolved stars; (iv) 
as a sink of heavy elements which if otherwise in the gas phase, would pro- 
foundly affect the interstellar gas chemistry; (v) as an efficient catalyst for 
the formation of H2 and other simple molecules as well as complex organic 
molecules (and as a protector by shielding them from photodissociating ul- 
traviolet [UV] photons) in the interstellar medium (ISM); (vi) as an efficient 
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agent for heating the interstellar gas by providing photoclectrons; (vii) as an 
important coolant in dense regions by radiating infrared (IR) photons (which 
is particularly important for the process of star formation in dense clouds 
by removing the gravitational energy of collapsing clouds and allowing star 
formation to take place); (viii) as an active participant in interstellar gas dy- 
namics by communicating radiation pressure from starlight to the gas, and 
providing coupling of the magnetic field to the gas in regions of low fractional 
ionization; (ix) as a building block in the formation of stars and planetary 
bodies and finally (x) as a diagnosis of the physical conditions (e.g. gas den- 
sity, temperature, radiation intensity, electron density, magnetic field) of the 
regions where dust is seen. 

The dust in the space between stars - interstellar dust - is the most 
extensively studied cosmic dust type, with circumstellar dust ("stardust"), 
comctary dust and IDPs coming second. Interstellar dust is an important 
constituent of the Milky Way and external galaxies. The presence of dust in 
galaxies limits our ability to interpret the local and distant Universe because 
dust extinction dims and reddens the galaxy light in the UV-optical-near-IR 
windows, where the vast majority of the astronomical data have been ob- 
tained. In order to infer the stellar content of a galaxy, or the history of star 
formation in the Universe, it is essential to correct for the effects of interstel- 
lar extinction. Dust absorbs starlight and reradiates at longer wavelengths. 
Nearly half of the bolometric luminosity of the local Universe is reprocessed 
by dust into the mid- and far-IR. 

Stardust, cometary dust and IDPs are directly or indirectly related to 
interstellar dust: stardust, condensed in the cool atmospheres of evolved stars 
or supernova ejecta and subsequently injected into the ISM, is considered as a 
major source of interstellar dust, although the bulk of interstellar dust is not 
really stardust but must have recondensed in the ISM pi]. Comets, because 
of their cold formation and cold storage, are considered as the most primitive 
objects in the solar system and best preserve the composition of the presolar 
molecular cloud among all solar system bodies. Greenberg [5D] argued that 
comets are made of unaltered pristine interstellar materials with only the most 
volatile components partially evaporated, although it has also been proposed 
that cometary materials have been subjected to evaporation, recondensation 
and other reprocessing in the protosolar nebula and therefore have lost all 
the records of the presolar molecular cloud out of which they have formed. 
Genuine presolar grains have been identified in IDPs and primitive meteorites 
based on their isotopic anomalies [8], indicating that stardust can survive 
journeys from its birth in stellar outflows and supernova explosions, through 
the ISM, the formation of the solar system, and its ultimate incorporation 
into asteroids and comets. 

Except in a few cases cosmic dust can be studied in situ (e.g. cometary 
dust [32], dust in the local interstellar cloud entering our solar system [22] ) 
or in terrestrial laboratories (e.g. IDPs and meteorites [8], cometary dust [7]), 
our knowledge about cosmic dust is mainly derived from its interaction with 
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electromagnetic radiation: extinction (scattering, absorption), polarization, 
and emission. Dust reveals its presence and physical and chemical properties 
and provides clues about the environment where it is found by scattering 
and absorbing starlight (or photons from other objects) and rcradiating the 
absorbed energy at longer wavelengths. 

This chapter deals with the optical properties of dust, i.e., how light is 
absorbed, scattered, and reradiated by cosmic dust. The subject of light scat- 
tering by small particles is a vast, fast-developing field. In this chapter I restrict 
myself to astrophysically-relevant topics. In 1 present a brief summary of 
the underlying physics of light scattering. The basic scattering terms are de- 
fined in Sj3j In fj4]I discuss the physical basis of the dielectric functions of dust 
materials. In Sj5]and fj6]I respectively summarize the analytic and numerical 
solutions for calculating the absorption and scattering parameters of dust. 

2 Scattering of Light by Dust: A Conceptual Overview 

When a dust grain, composed of discrete electric charges, is illuminated by 
an electromagnetic wave, the electric field of the incident electromagnetic 
wave will set the electric charges in the dust into oscillatory motion. These 
accelerated electric charges radiate electromagnetic energy in all directions, 
at the same frequency as that of the incident wave. This process, known 
as "scattering" (to be more precise, elastic scattering), removes energy from 
the incident beam of electromagnetic radiation. Absorption also arises as the 
excited charges transform part of the incident electromagnetic energy into 
thermal energy The combined effect of absorption and scattering, known as 
"extinction" , is the total energy loss of the incident wave [5] . 

The scattering of light by dust depends on the size, shape, and chemical 
composition of the dust and the direction at which the light is scattered. This 
can be qualitatively understood, as schematically shown in Figure [IJ by con- 
ceptually subdividing the dust into many small regions, in each of which a 
dipole moment will be induced when illuminated by an incident electromag- 
netic wave. These dipoles oscillate at the frequency of the incident wave and 
therefore scatter secondary radiation in all directions. The total scattered field 
of a given direction (e.g. Pi, P2) is the sum of the scattered wavelets, with 
their phase differences taken into account. Since these phase relations change 
for a different scattering direction, the scattered field varies with scattering 
direction^ The phase relations among the scattered wavelets also change with 
the size and shape of the dust. On the other hand, the amplitude and phase 
of the induced dipole moment for a given frequency depend on the material 
of which the dust is composed. Therefore, the scattered field is sensitive to 
both the size, shape and chemical composition of the dust [6]. 

1 An exception to this is the dust in the Rayleigh regime (i.e. with its size being 
much smaller than the wavelength) which scatters light nearly isotropically, with 
little variation with direction since for dust so small all the secondary wavelets 
are approximately in phase. 
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Fig. 1. A conceptual illustration of the scattering of light by dust. The dust is 
conceptually subdivided into many small regions. Upon illuminated by an incident 
electromagnetic wave, in each small region a dipole moment (of which the amplitude 
and phase depend on the composition the dust) will be induced. These dipoles 
oscillate at the frequency of the incident wave and scatter secondary radiation in all 
directions. The total scattered field of a given direction (e.g. Pi, P2) is the sum of the 
scattered wavelets, where due account is taken of their phase differences. Since the 
phase relations among the scattered wavelets change with the scattering direction 
and the size and shape of the dust, the scattering of light by dust depends on the 
size, shape, and chemical composition of the dust and the direction at which the 
light is scattered. 

3 Scattering of Light by Dust: Definitions 

As discussed in Sj2] when light impinges on a grain it is either scattered or 
absorbed. Let / (A) be the intensity of the incident light at wavelength A, the 
intensity of light scattered into a direction defined by 9 and <fi (see Fig. [2]) is 

m = «™ (!) 

where 0° < 9 < 180° is the scattering angle (the angle from the incident di- 
rection) , 0° < 4> < 360° is the azimuthal angle which uniquely determines the 
scattering plane defined by the incident direction and the scattering direc- 
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Fig. 2. Schematic scattering geometry of a dust grain in an incident radiation field 
of intensity J which scatters radiation of intensity 1(9, <j)) into a scattering angle 
9 (9 = 0°: forward scattering; 9 = 180°: backward scattering), an azimuthal angle 
4>, and a distance r from the dust. In a Cartesian coordinate system, the incident 
direction defines the +z axis. The scattering direction and the incident direction 
define the scattering plane. In the far-field region (i.e. fcr>l), I = I a F(9, (f>) /k 2 r 2 , 
where k = 2n/X is the wave number in vacuum. 



tion (sec Fig-EJ) HI F(@, 4 1 ) is the (dimensionlcss) angular scattering function, 
r 3> X/2n is the distance from the scatterer, and k = 2it/\ is the wave number 
in vacuum. The scattering cross section C sca , defined as the area on which the 
incident wave falls with the same amount of energy as that scattered in all 
directions by the dust, may be obtained by integrating the angular scattering 
distribution F(9,(f))/k 2 over all solid angles 

C S ca(A) = -j^J J F (9, (f>) sine dOdcj) , (2) 

2 When the scattering is along the incident direction (9 = 0°, i.e. "forward scat- 
tering") or the scattering is on the opposite direction of the incident direction 
(0 = 180°, i.e. "backward scattering), any plane containing the z axis is a suit- 
able scattering plane. 
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where F(9,<j))/k 2 (with a dimension of area)@ after normalized by C sca , is 
known as the phase function or scattering diagram 

_ F(g,0)/fc 2 
P{o,9) = 75 ■ (3) 

^sca 

The asymmetry parameter (or asymmetry factor) g is defined as the average 
cosine of the scattering angle 9 

g = (cos 9) = [ p(0,4>) cos9dH = -J- [ [ F(9, <t>) cos 9 sin 9 d9 d<t> . 

Ji-n K ^sca JO Jo 

(4) 

The asymmetry parameter 5, specifying the degree of scattering in the forward 
direction (9 = 0°), varies from —1 (i.e. all radiation is backward scattered like 
a "mirror") to 1 (for pure forward scattering). If a grain scatters more light 
toward the forward direction, g > 0; g < if the scattering is directed more 
toward the back direction; 5 = if it scatters light isotropically (e.g. small 
grains in the Raylcigh regime) or if the scattering is symmetric with respect 
to 9 = 90° (i.e. the scattered radiation is azimuthal independent and symmetric 
with respect to the plane perpendicular to the incident radiation). 

In radiative transfer modeling of dusty regions, astronomers often use the 
empirical Henyey-Greenstein phase function to represent the anisotropic scat- 
tering properties of dust [25] 

H (0) = J" 575 ' ( 5 ) 

4^(l + 5 2 -25COS#) 3/2 

Draine (2003) proposed a more general analytic form for the phase function 

H (9\ = - 3 ( 1 ~ g2 ) l + 77cos 2 

V[ ' ~ 4 7 r3 + 7? (l + 2.g 2 ) (i + g 2 2 g cos 9) 3/2 1 ' 

where 77 is an adjustable parameter. For 77 = Eq[6] reduces to the Henyey- 
Greenstein phase function. For 5 = and r\ = 1 EqjS] reduces to the phase 
function for Rayleigh scattering [T2"] . 

As discussed in £J5J both scattering and absorption (the sum of which is 
called extinction) remove energy from the incident beam. The extinction cross 
section, defined as 

total energy scattered and absorbed per unit time 

^ext — ^sca ~t~ ^abs '. T"j ] 7, . , , . j J 

incident energy per unit area per unit time 

is determined from the optical theorem which relates C ex t to the real part of 
the complex scattering amplitude S(9, (/>jf| in the forward direction alone [2"5] 



3 Also called the "differential scattering cross section" dC BC &/d£2 = F(6, 4>)/k 2 , it 
specifies the angular distribution of the scattered light [i.e. the amount of light (for 
unit incident irradiance) scattered into a unit solid angle about a given direction] . 

4 The angular scattering function F(8, (f>) is just the absolute square of the complex 
scattering amplitude S{0,(p): F(0,4>) = \S(9, 4>)\ 2 . 
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C cxt = ~—R c {S(0 = O°)} . (8) 

The absorption cross section C a bs is the area on which the incident wave falls 
with the same amount of energy as that absorbed inside the dust; C ex t, having 
a dimension of area, is the "effective" blocking area to the incident radiation 
(for grains much larger than the wavelength of the incident radiation, C ex t is 
about twice the geometrical blocking area) . For a grain (of size a and complex 
index of refraction m) in the Rayleigh limit (i.e. 2na/\<^i 1, 27ra|m|/A <C 1), 
the absorption cross section C a bs is much larger than the scattering cross 
section C sca and therefore C ox t ~ Cabs- Non-absorbing dust has C ex t = C sca . 

The albedo of a grain is defined as a = C sca /C ex t. For grains in the Rayleigh 
limit, a « since C sca "C C a b s - For Non-absorbing dust, a = 1. 

In addition to energy, light carries momentum of which the direction is that 
of propagation and the amount is hv/c (where h is the Planck constant, c is the 
speed of light, and v is the frequency of the light). Therefore, upon illuminated 
by an incident beam of light, dust will acquire momentum and a force called 
radiation pressure will be exerted on it in the direction of propagation of the 
incident light. The radiation pressure force is proportional to the net loss of 
the forward component of the momentum of the incident beam. While the 
momentum of the absorbed light (which is in the forward direction) will all 
be transfered to the dust, the forward component of the momentum of the 
scattered light will not be removed from the incident beam. Therefore, the 
radiation pressure force exerted on the dust is 

F pr = I C pr /c , C pr = C abs + (1 - g) C sca , (9) 

where C pr is the radiation pressure cross section, and I is the intensity (ir- 
radiance) of the incident light. 

In literature, one often encounters Q ex t, Qsca,, Q a bs, and Q w - the ex- 
tinction, scattering, absorption and radiation pressure efficiencies. They are 
defined as the extinction, scattering, absorption, and radiation pressure cross 
sections divided by the geometrical cross section of the dust C geo , 

n _ C ext _ _ C sca _ C a b s n _ C pr 

^gco ^gco ^gco ^gco 

For spherical grains of radii a, C goo ~na 2 . For non-spherical grains, there is 
no uniformity in choosing C goo . A reasonable choice is the geometrical cross 

section of an "equal volume sphere" C ge o = tt (3V/4ir) 2 ^ 3 f=a 1.21 F 2 / 3 where 
V is the volume of the non-spherical dust. 



4 Scattering of Light by Dust: Dielectric Functions 

The light scattering properties of dust are usually evaluated based on its opti- 
cal properties (i.e. dielectric functions or indices of refraction) and geometry 
(i.e. size and shape) by solving the Maxwell equations 
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VxE + --=0;VxH — = —J ; V • D = Anp ; V • B = ; (11) 

c at c at c 

where E is the electric field, B is the magnetic induction, H is the magnetic 
field, D is the electric displacement, J is the electric current density, and p 
is the electric charge density. They are supplemented with the constitutive 
relations (or "material relations") 

J = <tE ; D = eE = (1 + 4tt X ) E ; B = /iH ; (12) 

where a is the electric conductivity , s = 1 + 47r% is the dielectric function (or 
permittivity or dielectric permeability), % is the electric susceptibility, and 
/i is the magnetic permeability. For time-harmonic fields E, H oc cxp (— i tot), 
the Maxwell equations are reduced to the Helmholtz wave equations 

V 2 E + fc 2 E = ; V 2 H + fc 2 H = ; (13) 

where k — muj/c is the complex wavenumbcr, m = y//(e~-j-i Ana/uj) is the 
complex refractive index, and u> = 2nc/\ is the circular frequency. These equa- 
tions should be considered for the field outside the dust (which is the superpo- 
sition of the incident field and the scattered field) and the field inside the dust, 
together with the boundary conditions (i.e. any tangential and normal compo- 
nents of E are continuous across the dust boundary) . For non-magnetic dust 
(fi = 1), the complex refractive index is to = \J e + i 4-7T a/co. For non-magnetic, 
non-conducting dust, to = y/e. For highly-conducting dust, Aita /cu^$> 1, there- 
fore, both the real part and the imaginary part of the index of refraction are 
approximately \J2~ku jui at long wavelengths. 

The complex refractive indices to or dielectric functions e of dust are of- 
ten called optical constants, although they are not constant but vary with 
wavelengths. They are of great importance in studying the absorption and 
scattering of light by dust. They are written in the form of to = m! + im" 
and e = e' + ie". The sign of the imaginary part of m or e is opposite to 
that of the time-dependent term of the harmonically variable fields [i.e. to", 
e" >0 for E, Hex cxp (—icot); to", e" <0 for E, H oc exp (i u)t)]. The imagi- 
nary part of the index of refraction characterizes the attenuation of the wave 
(47rm"/A is called the absorption coefficient), while the real part determines 
the phase velocity (c/m 1 ) of the wave in the medium: for an electric field prop- 
agating in an absorbing medium of m — to' + im", say, in the x direction, 
E oc cxp(— Lum"x/c) exp[— iuj(t — m'x/c)]. 

The physical basis of the dielectric function e can readily be understood in 
terms of the classical Lorentz harmonic oscillator model (for insulators) and 
the Drude model (for free-electron metals). In the Lorentz oscillator model, 
the bound electrons and ions of a dust grain are treated as simple harmonic 
oscillators subject to the driving force of an applied electromagnetic field. The 
applied field distorts the charge distribution and therefore produces an induced 
dipole moment. To estimate the induced moments we consider a (polarizable) 
grain as a collection of identical, independent, isotropic, harmonic oscillators 
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with mass m and charge q; each oscillator is under the action of three forces: 
(i) a restoring force -ifx, where K is the force constant (i.e. stiffness) of the 
"spring" to which the bound charges are attached, and x is the displacement 
of the bound charges from their equilibrium; (ii) a damping force — 6x, where 
b is the damping constant; and (iii) a driving force gE produced by the (local) 
electric field E. The equation of motion of the oscillators is 

rax + bx + Kx = qF, . (14) 

For time harmonic electric fields E oc cxp (— iuit), we solve for the displacement 

(q/m) E 



j 2 — l^/LU 



o Q = K/m , 7 = b/m , (15) 



where tu is the frequency of oscillation about equilibrium (i.e. the restoring 
force is —moj 2 x). The induced dipole moment p of an oscillator is p = qx. Let 
n be the number of oscillators per unit volume. The polarization (i.e. dipole 
moment per unit volume) P = np = nqx is 

1 ,., 2 

E , w p 2 = Annq 2 /m , (16) 



4-7T 0J o 2 — OJ 2 — ijoj 



where lo v is the plasma frequency. Since P = %E = (e — 1) /4-7T E, the dielectric 
function for a one-oscillator model around a resonance frequency uj is 



£ = 1 + 4ttx = 1 + ^ \ — (17) 

, , ^P 2 (^o 2 — W 2 ) , , 

e' = 1 + 4^ x = 1 + p V 2 1— (18) 

(ui 2 — Ul 2 ) + 7 2 LJ 2 

e"=^ X " = • (19) 

(cjo — U ) + J Ld 

We see from EqfI2]e" > for all frequencies. This is true for materials close 
to thermodynamic equilibrium, except those with population inversions. If 
there are many oscillators of different frequencies, the dielectric function for 
a multiple-oscillator model is 

£ = 1+V 5 ^ . (20) 

where o> p j, u>j and jj arc respectively the plasma frequency, the resonant 
frequency, and the damping constant of the j-th oscillator. 

The optical properties associated with free electrons are described by the 
Drude model. The free electrons experience no restoring forces when driven by 
the electric field of a light wave and do not have natural resonant frequencies 
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Fig. 3. Refractive indices m' (upper panel), m" (lower panel) of neutral silicon 
nanoparticles (Si; solid lines), singly, positively charged (Z — +1) silicon nanopar- 
ticles of size a — 10 A (dashed lines), and SiC>2 glass (dot-dashed lines). Compared 
to neutral Si, in charged Si free electrons (if negative charged) or holes (if positively 
charged) contribute to the dielectric function. Crystalline Si is IR inactive since its 
lattice vibrations have no dipole moment; the bands at 6.91, 7.03, 7.68, 8.9, 11.2, 
13.5, 14.5, 16.4, and 17.9/im are due to multi-phonon processes. Taken from [48] . 



(i.e. u>o = 0). The Drudc model for metals is obtained directly from the Lorentz 
model for insulators simply by setting the restoring force in Eq |14l equal to 
zero. The dielectric function for free electrons is 

U) 2 

s = 1 , uj Pte 2 = 4irn e e 2 /m e (21) 

e' = l-^Z_ 2 (22) 

£ " = HTT^ (23) 

where Lo v>e is the plasma frequency, n e is the density of free electrons and 
m e is the effective mass of an electron. The damping constant 7 e = 1 /r e is 



Optical Properties of Dust 



11 



the reciprocal of the mean free time between collisions (r e ) which are often 
determined by electron-phonon scattering (i.e. interaction of the electrons 
with lattice vibrations) |f| For dust materials (e.g. graphite) containing both 
bound charges and free electrons, the dielectric function is 



2 ,,,2 . 

(24) 



The optical response of free electrons in metals can also be understood 
in terms of the electric current density J and conductivity <j{u>). The free 
electrons in metals move between molecules. In the absence of an external 
field, they move in a random manner and hence they do not give rise to a net 
current flow. When an external field is applied, the free electrons acquire an 
additional velocity and their motion becomes more orderly which gives rise 
to an induced current flow. The current density J = — n e ex is obtained by 
solving the equation of motion of the free electrons to c x = — eE — m e x/r e 

i= -r £ e 1 E 

m e 1 — i u)T e 

Since J = crE, the a.c. conductivity is a = a a j (1 — iujT e ), where a = 
n e T e e 2 1 'm e is the d.c. conductivity. The dielectric function for a free-electron 
metal is therefore 

e= l + i±g = l- a V / . (26) 

The real and imaginary parts of the dielectric function are not indepen- 
dent. They are related through the Kramers- Kronig (or dispersion) relations 

2„ f°° are" far) , ,„ . -2w „ f°° e'(x) , 
e'(u)=l + -P „ y ' dx , e"(w) = P , \ dx , (27) 

TT Jq X z — U! z TT Jq X Z — UJ" 

where P is the Cauchy Principal value of the integral 



5 For nano-sized metallic dust of size a (which is smaller than the mean free path of 
conduction electrons in the bulk metal), r c and 7 e are increased because of addi- 
tional collisions with the boundary of the dust: j e = jhuik+VF / where 7b u ik is 
the bulk metal damping constant, vf is the electron velocity at the Fermi surface, 
? is a dimensionless constant of order unity which depends on the character of the 
scattering at the boundary (?o is the effective mean free path for collisions with 
the boundary): ?=1 for classic isotropic scattering, ? = 4/3 for classic diffusive 
scattering, ? = 1.16 or ?=1.33 for scattering based on the quantum particle-in- 
a-box model (see [5] and references therein). Since oj 2 2>7^ in metals near the 
plasma frequency, e" can be written as e" = e"buik + kfw p 2 /(?ou 3 ). This, known 
as the "electron mean free path limitation" effect, indicates that for a metallic 
grain e" increases as the grain becomes smaller. 
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,r-^ d ,. ln [r-^ i+ r-^j. (28) 

J X z - UJ Z a^O IJ Q X Z - UJ Z J u+a X Z - UJ Z 

The real and imaginary parts of the index of refraction are also connected 
through the Kramers-Kronig relation. This also holds for the real and imagi- 
nary parts of the electric susceptibility. 

2„ f°° xm"(x) , ... , -2w„ f°° m'(x) , 
m'(w) = l + -P / 9 K t dx , m"{uj) = P / . y ' dx , (29) 

TT J Q X Z - U Z TT J Q X Z - U) Z 

x>) _ i P r «w 4 , x>) . ^ P r 4^ iix . (30) 

it Jo x ^ n Jo x — w 

The Kramers-Kronig relation can be used to relate the extinction cross 
section integrated over the entire wavelength range to the dust volume V 

/>OC 

/ C^ t {\)d\^iTx 2 FV , (31) 
Jo 

where F, a dimensionless factor, is the orientationally-averaged polarizability 
relative to the polarizability of an equal- volume conducting sphere, depending 
only upon the grain shape and the static (zero-frequency) dielectric constant 
e Q of the grain material [B3j . This has also been used to place constraints on 
interstellar grain models based on the interstellar depletions [46l [50] and the 
carrier of the mysterious 21 /xm emission feature seen in 12 protoplanetary 
nebulae [44] . 

For illustration, we show in Figure [3] the refractive indices of glassy SiC>2 
and neutral and singly-charged silicon nanoparticles (SNPs)0 The optical 
properties of SNPs depend on whether any free electrons or holes are present. 
The contribution of free electrons or holes to the dielectric function of SNPs 
is approximated by Se « — w p 2 r 2 / (w 2 r 2 + iuor), ui p 2 = 3\Z\ e 2 /a 3 m e ff, where 
e is the proton charge, Ze is the grain charge, a is the grain radius, r m cl/vf 
is the collision time (we take vf ~ 10 8 cms -1 ), and m e g is the effective mass 
of a free electron or hole. In Figure [3] the charged SNPs are taken to contain 
only one hole (i.e. Z = +i). We take m e g « 0.2 m e . 

5 Scattering of Light by Dust: Analytic Solutions 

For dust with sizes much smaller than the wavelength of the incident radia- 
tion, analytic solutions to the light scattering problem exist for certain shapes. 
Let a be the characteristic length of the dust, and x = 2Tra/X be the dimen- 
sionless size parameter. Under the condition of x^l and |m:r|<Cl (i.e. in 
the "Rayleigh" regime), C a b s = 47rfc Im{c*}, C sca = (87r/3) k A |a| 2 , where a. is 

6 SNPs were proposed as the carrier of the "extended red emission" (ERE), a broad, 
featureless emission band between ~ 5400 and 9000 A seen in a wide variety of 
dusty environments [721141] (but see [48] )• SNPs were also suggested to be present 
in the inner corona of the Sun [21] (but see [541 155] ) . 
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the complex electric polarizability of the dust. Apparently, C sca ^;C a bs and 
C cx t ~C a bs- In general, a is a diagonalized tensor Q for homogeneous spheres 
composed of an isotropic material, it is independent of direction 

We-1 

a " = ^—2 ' (32) 

where V is the dust volume^ For a homogeneous, isotropic ellipsoid, the 
polarizability for electric field vector parallel to its principal axis j is 

V e - 1 



4tt (e - l)L 3 + 1 



(35) 



where Lj is the "depolarization factor" along principal axis j (see PS])- The 
electric polarizability a is also known for concentric core-mantle spheres |69j , 
confocal core-mantle ellipsoids [131 Q35] , and multi-layered ellipsoids of equal- 
eccentricity [18j . For a thin conducting cylindrical rod with length 21 and 
radius r <C £, the polarizability along the axis of the rod is [38] 

I 3 

(36) 



31og(4Z/r )-7 



In astronomical modeling, the most commonly invoked grain shapes are 
spheres and spheroids (oblates or prolatcs)@ In the Raylcigh regime, their 
absorption and scattering properties arc readily obtained from Eqs. (|32I35[1 . 
For both dielectric and conducting spheres (as long as x<ti 1 and \mx\ <C 1) 



Cahs/V — X ^> C s , 

c (e' + 2) 2 + e" 2 



3 


© 4 


e-l 






e + 2 



(37) 



a can be diagonalized by appropriate choice of Cartesian coordinate system. It 
describes the linear response of a dust grain to applied electric field E: p=aE 
where p is the induced electric dipole moment. 

For a dielectric sphere with dielectric function given in Eq. ()17l) . in the Rayleigh 
regime the absorption cross section is 

C* hs (")/V=- J£l^L_ . (33) 

Similarly, for a metallic sphere with dielectric function given in Eg . (121 [I . 

It is seen that the frequency-dependent absorption cross section for both dielectric 
and metallic spheres is a Drude function. This is also true for ellipsoids. 
Spheroids are a special class of ellipsoids. Let r a , rj, and r c be the semi-axes of an 
ellipsoid. For spheroids, rt = r c . Prolates with r a >r\, are generated by rotating 
an ellipse (of semi-major axis r a and semi-minor axis r&) about its major axis; 
oblates with r a < rb are generated by rotating an ellipse (of semi-minor axis r a 
and semi-major axis rb) about its minor axis. 
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At long wavelengths, for dielectric dust e" oc to while e' approaches a con- 
stant much larger than e" (see Eqs ll8ll9")) . we see C a b s oc oje" oc oj 2 ; for 
metallic dust, e" oc \ju) while e' approaches a constant much smaller than e" 
(see Eqs l22l23"]) . we see C a b s oc u/e" oc w 2 ; therefore, for both dielectric and 
metallic dust C a b s oc A~ 2 at long wavelengths f*°l 

It is also seen from Eq. (f3"T)) that for spherical dust in the Rayleigh regime 
the albedo a f=a 0, and the radiation cross section C pr ks C a bs- This has an 
interesting implication. Let /3 pr (a) be the ratio of the radiation pressure force 
to the gravity of a spherical grain of radius a in the solar system or in debris 
disks illuminated by stars (of radius R+ and mass M*) with a stellar flux of 
F£ at the top of the atmosphere [53] , 

a ( x 3 E 2 / F* [q bs (a, A) + (1 - g) C sca («, A)] riA 

PprW = 77 3 (.OB J 

where G is the gravitational constant, and pdust is the mass density of the 
dust. For grains in the Rayleigh regime (<7« 0; C sca <C C a b s ; C abs oc a 3 ), we see 
/3 pr oc Cabs/a 3 is independent of the grain size a (see Fig.|4}! 

Spheroids are often invoked to model the interstellar polarization. In the 
Rayleigh approximation, their absorption cross sections for li ght polarized 
parallel (||) or perpendicular (_L) to the grain symmetry axis arq 11 ! 

,||, J- <^ 



g r/y=-lm( (e _ 1)J . ^ 09) 



where the depolarization factors parallel {L\\) or perpendicular (L±) to the 
grain symmetry axis are not independent, but related to each other through 
in + 2L ± = 1, with 



4m(l±|)-i] , e e = Vi-(W, a ) 2 (40) 



However, various astronomical data suggest a flatter wavelength-dependence (i.e. 
C a bs oc A - * 3 with j3<2): /?<2 in the far-IR/submm wavelength range has been 
reported for interstellar molecular clouds, circumstellar disks around young stars, 
and circumstellar envelopes around evolved stars. Laboratory measurements have 
also found (3 < 2 for certain cosmic dust analogues. In literature, the flatter (/3<2) 
long-wavelength opacity law is commonly attributed to grain growth by coagu- 
lation of small dust into large fluffy aggregates (see [47] and references therein). 
However, as shown in Eq. (|3ip . the Kramers-Kronig relation requires that f3 should 
be larger than 1 for A — +oo since F is a finite number and the integration in the 
left-hand-side of Eq. (|31|l should be convergent although we cannot rule out f3< 1 
over a finite range of wavelengths. 

For grains spinning around the principal axis of the largest moment of inertia, the 
polarization cross sections are C po i = (C^ bs — C a l bs )/2 for prolates, and C po i = 
(Cabs ~ C a 'bs) ^ or °blates; the absorption cross sections for randomly-oriented 
spheroids are C* a bs = (Cf bs + 2C a L bs )/3 05]. 
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Fig. 4. /3 P r - Ratio of radiation pressure to gravity for compact silicate grains in 
the debris disk around the Sun-like star BD+20 307 (GO, age ~ 300 Myr) . We note 
for nano-sized grains /3 pr ~0.12, independent of grain size; for grains larger than 
~ 0.3 /mi, /3 pr is inverse proportional to grain size. Song et al. [68] attributed the 
dust in this disk to recent extreme collisions between asteroids. Taken from |51| . 

for prolatcs (r a > rt,) where £ e is the eccentricity, and 



for oblates (r a <ri,). For spheres £|| = L± = l/3 and £ e = 0. For extremely 

elongated prolates or "needles" (r a ^> r^), it is apparent Cj^ s <C c| bs , we thus 
obtain 



where Lu ps (rb/r a ) 2 \n{r a /r^)- For dielectric needles, C a bs oc uje" oc A~ 2 at 
long wavelengths since L\\(e' - 1) + 1 > ^|| e " ( see S3); f° r metallic needles, 
for a given value of e" one can always find a sufficiently long needle with 
L\\e" < 1 and L\\{e' - 1) < 1 so that C a bs oc toe" oc a which can be very 
large (see [H]). Because of their unique optical properties, metallic needles 
with high electrical conductivities (e.g. iron needles, graphite whiskers) are 
resorted to explain a wide variety of astrophysical phenomena: (1) as a source 
of starlight opacity to create a non-cosmological microwave background by the 
thermalization of starlight in a steady-state cosmology [27]; (2) as a source of 




(41) 




(42) 
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the grey opacity needed to explain the observed redshift-magnitude relation of 
Type la supernovae without invoking a positive cosmological constant [T] ; (3) 
as the source for the submm excess observed in the Cas A supernova remnant 
[16] ; and (4) as an explanation for the flat 3-8 /im extinction observed for lines 
of sight toward the Galactic Center and in the Galactic plane [15]. However, 
caution should be taken in using Eq l42l (i.e. the Rayleigh approximation) 
since the Rayleigh criterion 2nr a \m\/ \ <C 1 is often not satisfied for highly 
conducting needles (see [45])F^I 

In astronomical spectroscopy modeling, the continuous distribution of el- 
lipsoid (CDE) shapes has been widely used to approximate the spectra of 
irregular dust grains by averaging over all ellipsoidal shape parameters [6]. 
In the Rayleigh limit, this approach, assuming that all ellipsoidal shapes are 
equally probable, has a simple expression for the average cross section 

(C ahs /V) = ^Im j-^-Logej , Loge = In \J s' 2 + e" 2 + i arctan (e"/e') 

(43) 

where Loge is the principal value of the logarithm of e. The CDE approach, 
resulting in a significantly-broadened spectral band (but with its maximum re- 
duced), seems to fit the experimental absorption spectra of solids better than 
Mie theory Although the CDE may indeed represent a distribution of shape 
factors caused either by highly irregular dust shapes or by clustering of spher- 
ical grains into irregular agglomerates, one should caution that the shape dis- 
tribution of cosmic dust does not seem likely to resemble the CDE, which as- 
sumes that extreme shapes like needles and disks are equally probable. A more 
reasonable shape distribution function would be like dP/dL\\ = 12L|| (1 — £||) 2 
which peaks at spheres (L" = 1/3). This function is symmetric about spheres 
with respect to eccentricity e and drops to zero for the extreme cases: infinitely 
thin needles (e — > 1, Lu^> 0) or infinitely flattened pancakes (e — * oo, L\\ — » 1). 
Averaging over the shape distribution, the resultant absorption cross section is 
Cabs = Jq dL^dP/dL^C a bs(L\\) where C a b s (£||) is the absorption cross section 
of a particular shape L\\ [SU [HH 05] • Alternatively, Fabian et al. [T7] proposed 
a quadratic weighting for the shape distribution, "with near-spherical shapes 
being most probable" . 

When a dust grain is very large compared with the wavelength, the electro- 
magnetic radiation may be treated by geometric optics: Q C xt = C e xt/C gco — > 2 



The "antenna theory" has been applied for conducting needle- like dust to estimate 
its absorption cross sections [75]. Let it be represented by a circular cylinder of 
radius r a and length / (r„<I). Let pr be its resistivity. The absorption cross 
section is given by C a bs = (47r/3c)(7rr^i/pR), with a long wavelength cutoff of A D = 
pnc(l/r a ) /ln(i/r a ) , and a short-wavelength cutoff of Amin ~ (2Trcm e ) / (pnn e e ), 
where m e , e, and n e are respectively the mass, charge, and number density of the 
charge-carrying electrons. 
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if x = 2-Ka/\^> \ and \m — \\x^$> For these grains C a b s ~ C geo ), 

the ratio of the radiation pressure to gravity (3 pr oc C a b s /o 3 oc 1/a. This is 
demonstrated in Figured] For dust with x^>l and \m — l|a:<Cl, one can 
use the "anomalous diffraction" theory [B5]. For dust with \m — l|a:<Cl and 
\m — 1| <C 1, one can use the Rayleigh-Gans approximation^ to obtain the 
absorption and scattering cross sections (6] [39] [69] : 

8 t r i 32\m-l\ 2 x A fAA . 

Q ab8 « -lm{mx} , Q sca » 27 - mx2 ■ (44) 

It is important to note that the Rayleigh-Gans approximation is invalid for 
modeling the X-ray scattering by interstellar dust at energies below 1 kcV. 
This approximation systematically and substantially overestimates the inten- 
sity of the X-ray halo below 1 keV [66] ■ 

6 Scattering of Light by Dust: Numerical Techniques 

While simple analytic expressions exist for the scattering and absorption prop- 
erties of dust grains which are either very small or very large compared to the 
wavelength of the incident radiation (see ^3]), however, in many astrophysi- 
cal applications we are concerned with grains which arc neither very small 
nor very large compared to the wavelength. Moreover, cosmic dust would in 
general be expected to have non-spherical, irregular shapes. 

Our ability to compute scattering and absorption cross sections for non- 
sphcrical particles is extremely limited. So far, exact solutions of scattering 



At a first glance, Q cx t — > 2 appears to contradict "common sense" by implying 
that a large grain removes twice the energy that is incident on it! This actually 
can be readily understood in terms of basic optics principles: (1) on one hand, all 
rays impinging on the dust are either scattered or absorbed. This gives rise to a 
contribution of C seo to the extinction cross section. (2) On the other hand, all the 
rays in the field which do not hit the dust give rise to a diffraction pattern that 
is, by Babinet's principle, identical to the diffraction through a hole of area C geo - 
If the detection excludes this diffracted light then an additional contribution of 
Cgeo is made to the total extinction cross section [6]. 

The conditions for the Rayleigh-Gans approximation to be valid are |m — 1| <C 1 
and \m — l|a:<Cl- The former ensures that the reflection from the surface of 
the dust is negligible (i.e. the impinging light enters the dust instead of being 
reflected); the latter ensures that the phase of the incident wave is not shifted 
inside the dust. For sufficiently small scattering angles, it is therefore possible for 
the waves scattered throughout the dust to add coherently. The intensity (I) of 
the scattered waves is proportional to the number (N) of scattering sites squared: 
I oc N 2 rxp 2 a 6 (where p is the mass density of the dust). This is why the X-ray 
halos (usually within ~ 1° surrounding a distant X-ray point source; [B3]) created 
by the small-angle scattering of X-rays by interstellar dust are often used to probe 
the size (particularly the large size end; [141 1671 [73] ). morphology (compact or 
porous; [581 166] ). composition, and spatial distribution of dust. 
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problems exist only for bare or layered spherical grains ("Mie theory"; [5]), 
infinite cylinders [53], and spheroids [3J [H [70]. The "T-matrix" (transition 
matrix) method, originally developed by Barber & Yeh [5] and substantially 
extended by Mishchcnko, Travis, & Mackowski EH], is able to treat axisym- 
metric (spheroidal or finite cylindrical) grains with sizes comparable to the 
wavelength. The discrete dipole approximation (DDA), originally developed 
by Purcell & Pennypacker |65j and greatly improved by Draine [lOj , is a pow- 
erful technique for irregular heterogeneous grains with sizes as large as several 
times the wavelength. The VIEF (volume integration of electric fields) method 
developed by Hage & Greenberg [25], based on an integral representation of 
Maxwell's equations, is physically similar to the DDA method. The microwave 
analog methods originally developed by Greenberg, Pedersen & Pedersen [21] 
provide an effective experimental approach to complex particles [23j . 

Although interstellar grains arc obviously non-spherical as evidenced by 
the observed polarization of starlight, the assumption of spherical shapes (to- 
gether with the Bruggeman or the Maxwell-Garnctt effective medium theories 
for inhomogeneous grains; [6]) is usually sufficient in modeling the interstellar 
absorption, scattering and 1R (continuum) emission. For 1R polarization mod- 
eling, the dipole approximation for spheroidal grains is proven to be successful 
in many cases. 

The DDA method is highly recommended for studies of inhomogeneous 
grains and irregular grains such as cometary, interplanetary, and protoplane- 
tary dust particles which are expected to have a porous aggregate structure. 
Extensive investigations using the DDA method have been performed for the 
scattering, absorption, thermal 1R emission, and radiation pressure properties 
of fluffy aggregated dust (e.g. see[2[2![2S[30l[311[33[32[^ 
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